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The paper describes the recent development in the simulation of the ion-cyclotron 
acceleration of the plasma in the VASIMR experiment. The modeling is done using an 
improved EMIR code for RF field calculation together with particle trajectory code for 
plasma transport calculation. The simulation results correlate with experimental data on 
the plasma loading and predict higher ICRH performance for a higher density plasma 
target. These simulations assist in optimizing the ICRF antenna so as to achieve higher 
VASIMR efficiency. 
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Nomenclature 

= magnetic vector potential (Wb / m) 

= magnetic induction ( 0-1 Tesla) 

= speed of light ( 3x1 (f m / s) 

= electric field (Volt / m) 

= identity tensor 
= imaginary unit 
= current (Ampere) 

= current density (0 - 10 s Ampere / m 2 ) 

= dialectic tensor 
= wave number (m 1 ) 

= azimuthal mode number 
= particle mass (3.35xl0~ 27 kg for Deuterium ion) 

= plasma particle density (0 - 1&° m' 3 ) 

= power ( 24,000 Watt) 

= power density (W / m 2 ) 

= particle charge (1.6x 10' 19 Coulomb) 

= resistance (0.2 Ohm - for circuit) 

= radial coordinates (m) 

= temperature (0- 100 eV) 

= time (s) 

= flow velocity (16, 000 m/s for inlet) 

= ion velocity (m / s) 

= ion energy (0-100 eV) 

= computational particle weght 
= thermal speed (m / s) 
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x j 

= 

cell position (m) 

X; 

= 

ion position (m) 

Z 

= 

axial coordinate (m) 

e 

= 

electric permittivity ( 8.85 10' 12 * F / m) 

<P 

= 

azimuthal coordinate (rad) 

n 

= 

efficiency 

p 

= 

magnetic permeability (7.25 Iff 6 Henry / m) 

Q 

= 

computational domain 

a 

= 

plasma conductivity tensor (S / m) 

(0 

= 

RF frequency (7.2 70 7 rad / s) 

¥ 

= 

magnetic flux (Tesla m 2 ) 

Subscripts: 


0 

= 

vacuum, inlet 

ANT 

= 

ICRF antenna 

c 

= 

circuit 

e 

ss 

electron 

1C 

= 

ion cyclotron 

i 

■ss 

ion 

j 

= 

cell number 

k 

= 

particle number 

l 

= 

particle 

P 

= 

plasma 

1 

= 

orthogonal directions to vacuum magnetic induction B 0 in the ( r, n) plane 

// 

= 

parallel to the magnetic induction B 0 


Rh“ 


I. Introduction 

ecent experiments at Advanced Space Propulsion Laboratory (ASPL) of the NASA Johnson Space Center 


have demonstrated significant ion heating at the Ion Cyclotron Range of Frequencies (ICRF) in the Variable 
Specific Impulse Magnetoplasma Rocket (VASIMR) 1 device. The current VASIMR experiment (VX-25, Figure 1) 
consists of the plasma production section including a helicon antenna with up to 20 kW of RF power, and a plasma 
booster area with an ICRF antenna that has up to 1.5 kW of RF power, though a 10 kW ICRF upgrade is underway. 
The ICRF coupling efficiency is related to plasma resistance R p and circuit resistance R c by 


Tj, c = 100% R P /(R P + RJ. 


( 1 ) 


Thus, to improve the ICRF efficiency, the plasma resistance needs to be higher and the circuit resistance lower. 



Figure 1. The VX-25 system. 
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As predicted by theory and numerical simulations, better ICRF efficiency occurs when there is a higher density 
in the plasma target. Mathematical simulations of ICRF antenna performance are currently done by the EMIR 
code 2 , which assumes that the RF electric field and RF current density depend on all three space independent 
variables and time, while the plasma target characteristics (magnetic field, plasma density, velocity, temperature, 
energy and absorbed power density) are in steady state and depend only on radial and axial position. All variables 
are presented in this paper have SI units, except for temperature and energy, which are measured in electron-volts. 


II. Mathematical Model for RF-field 

In the EMIR code 2 , the RF electric field, E, RF magnetic field, B, and RF antenna current density, j ANT , are 
expanded in a Fourier series in azimuthal coordinate <j>. Harmonic dependence with respect to time t is also 
assumed, so that the fields and currents in cylindrical coordinates (r, <j>, z) can be expanded as, for example: 

( 2 ) 

Ml 


Here m is an azimuthal mode number and co is RF frequency. 

The RF fields are obtained by solving Maxwell’s equations, which, written in harmonic form for the electric 
field E, are 


-V x V x E + ( ar / c 2 )E + ( ko / c 2 e )j p = - ico/xj ^ , (3) 

where c is speed of light, eis electric permittivity, and ]X is magnetic permeability in vacuum. 

The plasma current density j P is related to the electric field by a collisional cold plasma conductivity tensor a : 
j p -a-E. Equation (3) can then be represented by a system of independent equations for E m as presented by Stix 3 : 


-e~ lm *V x V x E m e""* +( 0 ) 2 /c 2 )K ■ E m = -iO)flj „ , 


(4) 


where K = I +(i/m)6 is a plasma dielectric tensor: 




(X 


0 



0 



(5) 


Here, the system of coordinates (1, <p, ||), 1 and || denotes, respectively, the direction perpendicular and parallel to 
the vacuum magnetic field B 0 in the (r, z) plane. In Eq. (4), j m is the m-th mode of the current density applied by the 
antenna, having only azimuthal non-zero component. 

The plasma dielectric tensor K is chosen from two options: either a collisional cold plasma ( K e ), or a reduced 
order kinetic description (ft r ) suitable for parallel wave propagation 4 . In the cold plasma model 3 , the entries of the 
dielectric tensor depend on the plasma density n, on the vacuum magnetic field B 0 and on the driving frequency co. 

In the warm plasma model 4 , a kinetic dispersion relation for parallel propagation 5 is used: 

*, 2 C 2 =6> 2 + X ar pl (0A* ] ( *j, ) . (6) 

l=ej 


This is solved for the wave number ky, by considering the effect of the plasma flow velocity V and the electron and 
ion temperature 7} on the conductivity. Parallel propagation can be considered using the same form for the dielectric 
tensor as that used in the cold plasma model if nonlocal effects caused by perpendicular propagation can be 
neglected. This model also considers the Doppler shift for the case of collisionless plasma. A reduced order 
formulation 5 is used for deriving the dielectric tensor components as follows 


3 

American Institute of Aeronautics and Astronautics 



where is the slow wave root from Eq. (4). The A?( ) and B?( Jfcj, ) functions 4 * in Eqs. (6) and (7) depend on 

plasma density n, velocity V" and temperature T,. Also, w^t = ^2 q t T u /m, is thermal speed, and oj p j = qf n/(e nti) is 
the plasma frequency. 

In VX-25 simulations, reported in this paper, the electron part of the dielectric tensor is that of the cold plasma 
model, while the ion dielectric tensor comes from the warm plasma model. 


III. Calculation of Plasma Target 

The EMIR code uses the following input data: 1) vacuum magnetic field. Bo, 2) target plasma density, 3) RF 
antenna geometry, and 4) computation domain and grid description. The vacuum magnetic field is calculated by a 
finite difference method using the MagStat code 6 . The plasma density can be calculated analytically or numerically 
from the particle trajectory code ParTraj 8 . The plasma velocity can be calculated either analytically or from the 
particle trajectory code. Ion temperature can be assumed constant or calculated by the particle trajectory code. 
Electron temperature is assumed constant. Electron density and electron velocity are assumed to be the same as ion 
density and ion velocity, respectivly. 



Figure 2. ID axial magnetic induction, density and velocity profiles 
for the baseline target plasma. For z > the plasma density 

n(r, z) and velocity V\\(r,z) are calculated from the particle 
trajectory code. Before the antenna, EMIR code uses analytical 
expressions for the plasma density and velocity. 


A. Vacuum magnetic 
field/induction 

Calculation of the vacuum magnetic 
induction B 0 for the set of cylindrical 
coils by MagStat code was described 
elsewhere 6 . The MagStat code solves 
the following equation for the magnetic 
potential A = (0, A^ 0): 


d 1 d(n4„) 

dr r dr dr 2 


-Mrj 0 \ 


( 8 ) 


A finite difference algorithm with a non- 
uniform adaptive mesh and a fast 
iterative solver are used. From the 
magnetic potential, the magnetic 
induction B 0 is calculated as 


B o =(B r ,0,B l ) = 


„ 1 3 , . 

rA t ) 

dz r dr 

(9) 


A ID plot of B z (0,z) is shown in Fig. 2. This field is used both in RF code EMIR, as well as in particle code 
ParTraj. A 2D contour plot of the field B 0 is shown in Fig. 4 as well as in a few figures that follow. 


B. Analytical presentation for the plasma target 

left 

Ion density is calculated analytically in that part of the computational domain Q for which z < Z ANT . For z > 

Z^ant ’ i° n density is calculated by a particle trajectory code as described in the following section. 

An example of a ID plot for the density profile and parallel velocity profile used for the antenna design is shown 
in Fig. 2. The plasma was assumed to be composed of 100% ionized deuterium. These baseline target parameters 
were chosen using fits to data from a reciprocating Langmuir probe, the actual currents used in the magnet coils for 
the experiment, and the flow rate of the feed gas. Significant uncertainties existed in the experimental 


4 

American Institute of Aeronautics and Astronautics 




measurements of the electron temperature this time, so the flow velocity was chosen as a parameter for study within 


the reduced order model. 



Figure 3. Radial density profiles were considered 
based on probe measurements and mapping back 
along magnetic field lines. 


The radial density profile measured at the probe 
location for the baseline study is shown in Fig. 3. The 
two-dimensional fit of the baseline density using these 
parameters is shown in Fig. 4. The density is assumed 
to be proportional to the magnetic induction strength 
along the axis, for the case where the ICRF antenna 
does not accelerate plasma too much. 

The general two-dimensional dependence of 
plasma density n(r,z) is reduced to a. simpler 
dependence n,{ yr), where yr = B ( /r,z)f is the magnetic 
flux. There are two options available to define the 
radial profile of the plasma density. In the 
“PARABOLIC” option the density is calculated 
through the following formula 
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In the “GAUSSIAN” option the density is calculated through the following formula 


nQf)=n 0 exp 
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( 11 ) 


The ion velocity and temperature fields for z > zj r . are defined by a particle trajectory code as described in 

the following section. For z < z) e J r the parallel plasma velocity and temperatures are assumed constant along the 
magnetic field lines as shown in Fig. 2. 


C. Particle trajectory simulation of the plasma target 

When plasma is accelerated by the RF power, the density is not proportional to the vacuum magnetic field. It is 
simulated by a simple version of the particle method, called the “particle trajectory method”, where the plasma fluid 
characteristics (density, velocity, temperature) are assumed time-independent. Also, collisions are neglected. 

The ion motion satisfies the following equation of motion with respect to ion position and velocity vectors, #,• = 
(x n x$, x 7 ) and v ( = (v„ v$, vj, as functions of time t: 


m, ~ = q, (v, xB(x,,t) + E (x i , t)), ~ 
dt dt 


• = v, . 


( 12 ) 


Single particle trajectories are integrated from Eq. (12) with an adaptive time-scheme, which can quickly solve 
extensive particle simulations for systems of hundreds of thousands of particles in a reasonable time (10 min), and 
without the need for a powerful supercomputer. The particle calculation method is described in previous 
publications 7 ' 8,9,10 . During extensive particle simulation, one needs to define an initial ion velocity distribution, as 
will be discussed in the following section. 

Notice, that implementation of the formula (12) requires interpolation of the pre-calculated fields B and E from 
their grids to the arbitrary position. This is implemented efficiently using a grid-table method. 
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D. Particle to Grid Weighting 

Here we present our approach of calculating plasma fluid characteristics (density n, velocity V, temperature T ; ) 
based on the kinetic theory, from the calculated particle trajectories, defined by particle positions and velocities 
(a:, and vj. The discrete ion density n is defined constant at each finite difference cell Xj, using the formula 

n(Xj) = w^Count(x k e X j ) , (13) 

k 

where x k is a position of ^-particle, w is a particle weight calculated, such that it makes the grid density equal given 
value at given point n(X 0 )=n°: w = n°J y £ i Count(x k e X 0 ). 

Now, to describe the method of trajectories 11 . All trajectories start at the same axial position. The choice of that 
location is made from the following consideration. Only the ICRF area is considered with non-collisional behavior, 
so that location should be downstream of the flow from the magnet 3, where maximal magnetic strength occurs. 

Also, it has to be upstream of the flow from the Ion Cyclotron heating location. Thus it can be chosen near z l ^j . 

First, the non-uniform grid is generated defining the finite-size cells Xj. Assuming that the vacuum magnetic 
induction B 0 is precomputed on B-grid, as well as RF field E m on RF-grid, allows grid tables to be generated. These 
grid tables are the fast way to find a grid cell corresponding to the given position x h 

With a given distribution for the initial position and velocity vectors at a z near z^nt > a l ar g e number (order of 

10 s ) of ion trajectories are calculated. For simplicity it is assumed that the distribution function at this initial 
location is a product of radial distribution, velocity magnitude distribution and two velocity slope distributions. 
Those distributions are chosen such way to satisfy the experimental data, like in Fig. 3. Current simulation assume 
that a Maxwellian distribution of ion velocity occurs in the flow before it goes into the computational domain inlet 

at z- zf‘ NT • Each single trajectory is used to generate many particles distributed along it with equal time steps 
between them. 

Other fluid quantities (velocity V, ion current density j, temperature T and energy W) are calculated by a 
technique similar to that used to calculate the ion density in (9), as presented below: 


V(Xj) = { v k : x t e : x k e X } } /%Count{x k e X ,) , = g,w^{v k : x* e X .} , (14) 

k / k k 

T ll (X j ) = m i '£{(v u -V [{ f:x k zX j } /( q l '2 J Count{x k e Z ; }1, W^Xj) = : ** e Xj) ji 2q j '£ j Count{x k e Xj] 1. 



trajectory code and magnetic configuration of the VX-25 experiment, 
method. 


After these 

summations are 

performed, they are 
smoothed to reduce 
noise. While Fig. 2 
demonstrates a ID plot 
of the n(0, z). Fig. 4 
shows a 2D plot of 
n(r,z). 

After getting the 
plasma density 

calculated by the 
particle trajectory code, 
the solution is expanded 
to the rest of the 
computational domain 
using an analytical 
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Figure 5. 2D contour plot of the plasma azimuthal velocity V/r,z) 
calculated by the particle trajectory code and magnetic 
configuration of the VX-25 experiment. 


Plasma velocity is then presented in 
the magnetic coordinates: (JL, <j>, ||). As 
the numerical experiment shows, and 
theory predicts, the first component of 
the velocity vector is negligible: V± s 0, 
and can be assumed zero. The azimuthal 
component of the velocity vector is zero 
on the axis but grows toward the edge of 
the plasma especially in the area of the 
antenna, as shown in Fig. 5. This 
demonstrates the interaction with ion- 
cyclotron electromagnetic waves. 

The large azimuthal velocity on the 
edge of the plasma can be explained by 
the diamagnetic drift, presented by 
analytical formula 


V D (r,z) = -V(nT 1 )xB/nB 2 , 


(15) 


which gives the following azimuthal component: 

V/r,z) = -VJnT ± )/nB. 


(16) 



calculated by the particle trajectory code. 


The parallel velocity Vjj does not change 
across the plasma radius; an example is 
shown in Fig. 2. 

The ion temperature is also best 
represented in magnetic coordinates (1, 
<j>, ||). As numerical experiment shows 
and theory predicts, the first and second 
temperatures are about the same: T ± = 
These and Tj| don’t change much across 
the plasma radius. This is why the 
temperature can be presented by ID plot 
as shown in Fig. 6. 

Figure 6 is a ID plot of the average 
energies and temperatures of ions along 
the plasma cross-section. The total 
energy |W| is a sum of the parallel and 
orthogonal energies and Wx- 


E. Diamagnetic field calculation 

To calculate the plasma magnetic field, the plasma current has to be derived first. In an axisymmetric system 
with cylindrical symmetry, this current is purely azimuthal: 


jj r,z) = qrf r,zW t ( r,z)~ - q ^ (17) 

D 

and it has a direction opposite to the electro-magnet current density jo- The plasma current generates a plasma 
magnetic field that satisfies Ampere’s law: 


VxB p =/ij p . 


( 18 ) 
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0.5 0.6 0.7 0.8 0.9 1 1.1 

Figure 7, 2D contour plots of the diamagnetic current density 
calculated by the particle trajectory code and the magnetic field 
lines for the corresponding diamagnetic plasma field Bp. 


field, so the reduction of the total magnetic field is not significant. 


Figure 7 demonstrates the plasma 
diamagnetic current and the magnetic 
field calculated by last equation. The 
diamagnetic effect is essential only for 
high magnetic field in the area close to 
the thruster core and becomes negligible 
further away from it. The plasma 
magnetic field is calculated using the 
same solver, as used for the vacuum 
magnetic field calculation. The only 
difference in this calculation is the 
presence of a current density source j p . 
Calculation of B p is iterated with the 
calculation of plasma velocity and 
density. For the VX-25 experiment, the 
maximal diamagnetic field has a strength 
less than 1% of the vacuum magnetic 


IV. RF Power Absorption Simulation 

The RF power absorbed by the plasma is calculated by the following formula 


P = Re 


[Ej p 'd£2 


^ m l=ej 


| E m -jfdQ 


^ m l=e,i 


\E m -(d r E m ) d£2 = Jp( r,z)d£2 , (19) 


where p( r,z ) = Re E'" • (&, -E m ) 1 is a RF power density. 

^ m l=e,i L J 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Figure 8. 2D plot of the magnetic induction configuration for the VASIMR 
experiment with the ICRF power density p(r,z)- Dashed lines indicate 
resonance location. 

resonance, shown by the dashed line indicating B a (r,z)= B ICDop (r, z) = (a> - 
calculated by the EMIR code and Vj| is calculated by the particle trajectory code. 


The new model in EMIR 
allows the simulation of 
wave-plasma interactions 
with greater physical 
accuracy then was possible 
in the previous cold plasma 
model because the Doppler 
shift 12 effect for ions is now 
included. As shown in Figs. 
8 and 9, the power density 
distribution vs. longitudinal 
coordinate (and therefore vs. 
axial magnetic field) 
features a peak shifted 
downstream from the cold 
plasma ion cyclotron 
resonance due to the 
Doppler effect. The cold 
resonance corresponds to 
the magnetic field value 
B 0 (r,z) = B lc =a)m i /q i . 
Most of the power gets 
absorbed by the plasma at 
the Doppler-shifted 

■ L(r,z)V ]l (r,z})m i /q 1 . k\ 


is 
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Figure 9. ID plots of the magnetic induction configuration for the VASIMR 
experiment with the ICRF power absorption calculated by the EMIR code. 
Vertical dash lines indicate the ICRF antenna boundaries. 


cold ICRF resonance at the location where B 0 (r, z) = B IC Dap ( r , z) . 


Figure 9 shows some ID 
plots of the magnetic 
induction configuration for 
the VASIMR experiment 
VX-25 with the ICRF 
power absorption calculated 
by the EMIR code. Vertical 
dashed lines indicate the 
ICRF antenna boundaries. 
The absorbed power density 
has three peaks. The first 
peak indicates the cold IC 
resonance upstream the 
plasma flow, and the second 
peak indicates the cold IC 
resonance downstream the 
plasma flow. Both cold 
resonance peaks are located 
at B 0 (r, z) = B lc , and B ic is 
marked by the horizontal 
dashed line. Due to the 
Doppler shift most of the 
total power gets absorbed 
by the plasma at the third 
peak downstream from the 


The RF power absorption by the plasma for a given antenna current determines the “plasma loading resistance”, 
R p , a figure of merit for the antenna design. In a lumped circuit model, the resistance for each antenna segment can 
be defined as twice the power emitted by that segment divided by the square of the current in that segment. In order 
to efficiently couple RF power to the plasma, the loading resistance for the entire antenna must be substantially 
larger than the vacuum loading resistance, which is due to the finite resistance of the entire antenna driving circuit. 

An important goal for of the ICRF antenna design is to achieve high RF plasma resistance. Since the absorbed 


RF power is proportional to the plasma resistance as 


Pant ~~Ep {jnax(Re( J ANT ))J , the EMIR code allows 


calculating the resistance from the absorbed power for any given antenna current. The plasma resistance R p , 
calculated from the absorbed RF power as R p = 2P ANT (max(Re( J ^ )))' is compared with the circuit loading R c to 
analyze the fraction of the ICRF power going into plasma. Based on the calculated plasma resistance, the ICRF 
efficiency can be calculated based on the formula: J] ICRF = 100% R p j [R p + R c ) . For the present configuration of the 

VX-25 experiment, the circuit resistance is measured as R c = 0.2 Ohm. 

The EMIR code normally runs with the amplitude of the antenna current set equal to 1 Ampere. The plasma 
resistance can be calculated from the absorbed power P p : R p = 2 P p /(1A f. Absorbed power is a sum of the powers 
absorbed by electrons and ions. Objectives of ICRF antenna design are the increase the ICRF efficiency, and the 
increase in ion power fraction. 

The RF electric field in the EMIR code is linear with respect to the antenna current j ANT . After EMIR calculates 
the RF field, with antenna current amplitude set equal to 1 Amp, the solution can be scaled to match the ICRF power 
we wish to simulate: 


E, c ~ Ej 


After scaling, the RF field can be used for simulating the plasma target using the particle trajectory code. 


( 20 ) 
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V. Parameter Scan Studies 
Extensive parameter studies are 
being conducted at the ASPL to 
optimize ICRF heating performance 
both in the VX-25 experiment as well as 
in the EMIR simulations. Figure 10 
demonstrates the effect of ion flow rate 
on ICRF efficiency and ion power 
fraction, by EMIR simulation for a 
Deuterium plasma. The plasma density 
was fixed, while the input flow velocity 
was varied. As shown in Fig. 10, the 
ICRF efficiency goes up with the ion 
flow velocity, while the ion power 
fraction does not change much. 

Figure 1 1 compares the experimental 
data with simulation results for the 
effect of ion flow rate on plasma 
loading for a Helium plasma. The 
simulations were conducted for 
different ion densities as well as for different flow velocities. Both experimental and simulation results agree that 
plasma loading goes up with the ion flow rate. 


Figure 10. 
rate. 


ICRF performance dependence on the input ion flow 
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Figure 11. Input gas flow scan results for experiment (left) and numerical simulation (right) 


100 


Figure 12 demonstrates EMIR results on plasma loading and ion power fraction for different antenna geometries. 
A higher number of straps makes the loading higher. The triple strap antenna seems optimal, since the quadruple 
strap does not increase ion loading much. The quarter twist antenna seems have better loading than the half twist 
antenna or antenna with no twist. 
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Double Strap 
Half Twist Anlenna: 

Ion loading 70 mf) 
Electron loading 90 mf) 
ICRF efficiency 38% 

Ion power fraction 58% 

Figure 12. 
geometry 


Triple Strap 

Half Twist Anlenna 
Ion loading 138 mi) 
Electron loading 96 mf) 
ICRF efficiency 54% 

Ion power fraction 59% 


Quadruple Strap 
Half Twist Antenna 
Ion loading 149 mf) 
Electron loading 197 mfl 
ICRF efficiency 63% 

Ion powerfraction 43% 




Half Twist Antenna 
Unloading 138 mf) 
Electron loading 96 mD 
ICRF efficiency 54% 

Ion power traction 59% 


Quarter Twist Antenna: 
Ion loading 152 mf) 
Electron loading 168 mf) 
ICRF efficiency 62% 

Ion power fraction 47% 


Zero Twist Antenna 
Ion loading 134 mf) 
Electron loading 262 mf) 
ICRF efficiency 66% 

Ion power fraction 34% 


Analysis of the ICRF plasma loading and the ICRF efficiency as a function of the antenna 


Figure 13 shows experimental and simulation results for ICRF loading for different magnetic field slopes under 
the ICRF antenna. Three magnetic field profiles were compared: flat (zero), slightly negative, and slightly positive 
slopes with the cold ICRF resonance under the middle of the antenna. As both experiment and simulations agree, 
the best ICRF performance is provided for a magnetic field profile with a slightly positive slope under the antenna. 
This profile corresponds to the multiple cold ICRF resonance under the antenna. Due to the Doppler shift, all 
slopes had maximal RF absorption at the area downstream from the antenna. 


+VB VB~0 -VB 




0.9 0.95 1 1.05 1.1 1.15 1.2 

CO/CDci 

Figure 13. Analysis of the ICRF plasma resistance as a function of the magnetic field slope under antenna 
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VI. Conclusions 

Effective ion cyclotron wave coupling to the plasma has been observed in the VASIMR experiment. Recent 
plasma parameters considered with the latest ICRF antenna design give a plasma loading with R p » R c , that is, 
most of the RF transmitter power gets absorbed by the plasma, with more than half directly coupled to the ions. 
Future design improvements are directed to increase overall VASIMR efficiency by minimizing RF circuit losses, 
increasing plasma flux from the helicon source and optimizing the magnetic geometry. 

The EMIR simulation of RF wave propagation in the VASIMR plasma allows us to design the ICRF antenna 
with improved efficiency, as well as to optimize the magnetic field profile. Recent developments in the EMIR code 
have resulted in a code performance improvement that will enable further refinement in the physics of the model 
without sacrificing the fast turnaround of simulation results that is required in the present design effort to support the 
VASIMR experiment. 

Future plasma heating simulations in the VASIMR system will involve integration of the EMIR code and particle 
simulation code into “Numerical VASIMR” simulation software. The plan is to have a self-consistent modeling tool 
of all significant VASIMR areas: plasma production, plasma booster and exhaust. The next version of the EMIR 
code will also be parallelizable, and the next version of the particle trajectory code will involve Monte-Carlo 
collisions. 
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